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We study the properties of quantum stabilizer codes that 
embed a finite-dimensional protected code space in an infinite- 
dimensional Hilbert space. The stabilizer group of such a code 
is associated with a symplectically integral lattice in the phase 
space of 2N canonical variables. From the existence of sym- 
plectically integral lattices with suitable properties, we infer a 
lower bound on the quantum capacity of the Gaussian quan- 
tum channel that matches the one-shot coherent information 
optimized over Gaussian input states. 



I. INTRODUCTION 



A central problem in quantum information theory is 
to determine the quantum capacity of a noisy quantum 
channel the maximum rate at which coherent quantum 
information can be transmitted through the channel and 
recovered with arbitrarily good fidelity A general 

solution to the corresponding problem for classical noisy 
channels was found by Shannon in the pioneering paper 
that launched classical information theory With 
the development of the theory of quantum error correc- 
tion ||^|, considerable progress has been made toward 
characterizing the quantum channel capacity 0, but it 
remains less well understood than the classical capacity. 

The asymptotic coherent information has been shown 
to provide an upper bound on the capacity and 
a matching lower bound has been conjectured, but not 
proven JTo| . Unfortunately, the coherent information is 
not subadditive H], so that its asymptotic value is not 
easily computed. Therefore, it has been possible to verify 
the coherent information conjecture in just a few simple 
cases Q. 

One quantum channel of considerable intrinsic inter- 
est is the Gaussian quantum channel, which might also 
be simple enough to be analytically tractable, thus pro- 
viding a fertile testing ground for the general theory of 
quantum capacities. A simple analytic formula for the 
capacity of the Gaussian classical channel was found by 
Shannon |3||I| . The Gaussian quantum channel was stud- 
ied by Holevo and Werner Jl3|, who computed the one- 
shot coherent information for Gaussian input states, and 
derived an upper bound on the quantum capacity. 
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Lower bounds on the quantum capacity of the Gaus- 
sian quantum channel were established by Gottesman, 
Kitaev and Preskill They developed quantum 

error-correcting codes that protect a finite-dimensional 
subspace of an infinite-dimensional Hilbert space, and 
showed that these codes can be used to transmit high- 
fidelity quantum information at a nonzero asymptotic 
rate. In this paper, we continue the study of the Gaus- 
sian quantum channel begun in Our main result is 
that the coherent information computed by Holevo and 
Werner is in fact an achievable rate. This result lends 
nontrivial support to the coherent information conjec- 
ture. 

We define the Gaussian quantum channel and review 
the results of Holevo and Werner in Sec. II. In Sec. 
Ill, we describe the stabilizer codes for continuous quan- 
tum variables introduced in 14 1, which are based on 
the concept of a symplectically integral lattice embed- 
ded in phase space. In Sec. IV and V, we apply these 
codes to the Gaussian quantum channel, and calculate an 
achievable rate arising from lattices that realize efficient 
packings of spheres in high dimensions. This achievable 
rate matches the one-shot coherent information Iq of the 
channel in cases where 2 Iq is an integer. Rates achieved 
with concatenated coding are calculated in Sec. VI; these 
fall short of the coherent information but come close. In 
Sec. VII, we consider the Gaussian classical channel, and 
again find that concatenated codes achieve rates close to 
the capacity. Sec. VIII contains some concluding com- 
ments about the quantum capacity of the Gaussian quan- 
tum channel. 



II. THE GAUSSIAN QUANTUM CHANNEL 

The Gaussian quantum channel is a natural general- 
ization of the Gaussian classical channel. In the classical 
case, we consider a channel such that the input x and 
the output y are real numbers. The channel applies a 
displacement to the input by distance £, 



(1) 



where £ is a Gaussian random variable with mean zero 
and variance a 2 ; the probability distribution governing £ 
is 

1 



p(0 



(2) 



Similarly, acting on a quantum system described by 
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canonical variables q and p that satisfy the commutation 
relation [q,p] = ih, we may consider a quantum channel 
that applies a phase-space displacement described by the 
unitary operator 

D(a) = exp (aa) + a* a) , (3) 

where a is a complex number, [a,a^] = 1, and q, p can 
be expressed in terms of a and a' as 




This quantum channel is Gaussian if a is a complex Gaus- 
sian random variable with mean zero and variance a 2 . 
In that case, the channel is the superoperator (trace- 
preserving completely positive map) £ that acts on the 
density operator p according to 

p -> £( P ) = -L / d 2 a e-W'* 2 D(a)pD(ay . (5) 
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In other words, the position q and momentum p are dis- 
placed independently, 

q -> q + £,q , V^P + ip , (6) 

where £ q and £ p are real Gaussian random variables with 
mean zero and variance a 2 = Ha 2 . 

To define the capacity, we consider a channel's nth 
extension. In the classical case, a message is transmitted 
consisting of the n real variables 

x = (xi,x 2 , ■ ..,x n ) , (7) 
and the channel applies the displacement 

x^x + t, f= (6, 6, ■■■,£«) , (8) 

where the £j's are independent Gaussian random vari- 
ables, each with mean zero and variance a 2 . A code 
consists of a finite number m of n-component input sig- 
nals 

x (a) , a= 1,2,..., m (9) 

and a decoding function that maps output vectors to the 
index set {1,2,..., m}. We refer to n as the length of the 
code. 

If the input vectors were unrestricted, then for fixed 
a 2 we could easily construct a code with an arbitrarily 
large number of signals m and a decoding function that 
correctly identifies the index (a) of the input with ar- 
bitrarily small probability of error; even for n = 1 we 
merely choose the distance between signals to be large 



compared to a. To obtain an interesting notion of capac- 
ity, we impose a constraint on the average power of the 
signal, 

± <p, do) 

i 

for each a. We say that a rate R (in bits) is achievable 
with power constraint P if the there is a sequence of 
codes satisfying the constraint such that the /3th code in 
the sequence contains mp signals with length rip, where 

R = lim — log 2 mp , (11) 

/3^oo Tip 

and the probability of a decoding error vanishes in the 
limit /3 — * co. The capacity of the channel with power 
constraint P is the supremum of all achievable rates. 

The need for a constraint on the signal power to de- 
fine the capacity of the Gaussian classical channel can 
be understood on dimensional grounds. The classical ca- 
pacity (in bits) is a dimensionless function of the variance 
a 2 , but a 2 has dimensions. Another quantity with the 
dimensions of a 2 is needed to construct a dimensionless 
variable, and the power P fills this role. 

In contrast, no power constraint is needed to define 
the quantum capacity of the quantum channel. Rather, 
Planck's constant h enables us to define a dimensionless 
variance a 2 = a 2 /h, and the capacity is a function of 
this quantity. In the quantum code consists of 

an encoding superoperator that maps an m-dimensional 
Hilbert space H. m into the infinite-dimensional Hilbert 
space Ti^ N of N canonical quantum systems, and a de- 
coding superoperator that maps 7^ 8JV back to Ti. m . We 
say that the rate R (in qubits) is achievable if there is a 
sequence of codes such that 

R = lim — log 2 mp , (12) 

/3^oo iV p 

where arbitrary states in TL m can be recovered with a 
fidelity that approaches 1 as f3 — * oo. The quantum 
capacity Cq of the channel is defined as the supremum 
of all achievable rates. 

Holevo and Werner |D| studied a more general Gaus- 
sian channel that includes damping or amplification as 
well as displacement. However, we will confine our at- 
tention in this paper to channels that apply only dis- 
placements. Holevo and Werner derived a general upper 
bound on the quantum capacity by exploiting the prop- 
erties of the "diamond norm" (norm of complete bound- 
edness) of a superoperator. The diamond norm is defined 
as follows: First we define the trace norm of an operator 
X as 

pfllte = tn/xtX , (13) 
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which for a self-adjoint operator is just the sum of the 
absolute values of the eigenvalues. Then a norm of a 
superoperator £ can be defined as 



\\£ || so = sup 



HAIL. 



(14) 



The superoperator norm is not stable with respect to 
appending an ancillary system on which £ acts trivially. 
Thus we define the diamond norm of £ as 



|£||o = sup \\£ (g> I, 



11 1 1 SO , 



(15) 



where /„ denotes the n-dimensional identity operator. 
(This supremum is always attained for some n no larger 
than the dimension of the Hilbert space on which £ acts.) 
Holevo and Werner showed that the quantum capacity 
obeys the upper bound 



C Q {£) <log 2 ||£oT| 



(16) 



where T is the transpose operation defined with respect 
to some basis. In the case of the Gaussian quantum chan- 
nel, they evaluated this expression, obtaining 

C Q (a 2 )<log 2 (h/a 2 ) (17) 

for h/a 2 > 1, and C Q {a 2 ) = for h/a 2 < 1. 

Holevo and Werner p3| also computed the coherent in- 
formation of the Gaussian quantum channel for a Gaus- 
sian input state. To define the coherent information of 
the channel £ with input density operator p, one intro- 
duces a reference system R and a purification of p, a pure 
state 1$) such that 



tr fl (|$)<$|)=p 



(18) 



Then the coherent information Iq is 

I Q (£,p) = S(£(p))-S(£®I R (\<f>){$\)) , (19) 

where S denotes the Von Neumann entropy, 

S(p)=-ti(plog 2 p) . (20) 

It is conjectured |l^J^,^| that the quantum capacity is 
related to the coherent information by 



where 



C Q (£)= lim --C n {£) , 

n— >oo Tl 



C n (£)= S upI Q (£® n ,p) 
p 



(21) 



(22) 



Unlike the mutual information that defines the classical 
capacity, the coherent information is not subadditive in 



general, and therefore the quantum capacity need not 
coincide with the "one-shot" capacity C\. Holevo and 
Werner showed that for the Gaussian quantum channel, 
the supremum of Iq over Gaussian input states is 



('Q)max = ^g 2 (h/ea 2 ) 



(23) 



(where e = 2.71828..) for h/ea 2 > 1, and (/ Q ) max = 
for h/ea 2 < 1. According to the coherent-information 
conjecture, eq. should be an achievable rate. 



III. QUANTUM ERROR CORRECTING CODES 
FOR CONTINUOUS QUANTUM VARIABLES 



The lattice codes developed in |l4| are stabilizer codes 
|p^| , ^6| that embed a finite-dimensional code space in 
the infinite-dimensional Hilbert space of TV "oscilla- 
tors," a system described by 2N canonical variables 
qi, q%, • • • q_N,pi,P2, ■ ■ ■ ,pn- That is, the code space is the 
simultaneous eigenstate of 2A commuting unitary oper- 
ators, the generators of the code's stabilizer group. Each 
stabilizer generator is a Weyl operator, a displacement in 
the 2iV-dimensional phase space. 

Such displacements can be parametrized by 2N real 
numbers ot\, a.2, ■ ■ ■ , un, 01,02, ■ ■ ■ , @N, and expressed as 



U (a, (3) = exp 




(24) 



Two such operators obey the commutation relation 

U(a,/3)U(a',/3') = e^^^'P^Uia' , f3')U(a, /3) , (25) 



where 



u(a/3,a'p') = a -13' -a' ■ (i 



(26) 



is the symplectic form. Thus Weyl operators commute if 
and only if their symplectic form is an integer. 

The 2N generators of a stabilizer code are commuting 
Weyl operators 



U 



(a (a) ,/3 (a) ) , a = 1,2,. 



,27V 



(27) 



Thus the elements of the stabilizer group are in one-to- 
one correspondence with the points of a lattice £ gener- 
ated by the 2iV vectors u (a) = (a (a) , /3 (a) ). These vectors 
can be assembled into the generator matrix M of C given 
by 



M 



(28) 



3 



Then the requirement that the stabilizer generators com- 
mute, through eq. (p5|), becomes the condition that the 
antisymmetric matrix 

A = MuiM T (29) 

has integral entries, where M T denotes the transpose of 
M, to is the 2JV x 2N matrix 

"=(-l o) (30 » 

and In is the N x N identity matrix. If the generator 
matrix M of a lattice C has the property that A is an 
integral matrix, then we will say that the lattice C is 
symplectically integral. 

Encoded operations that preserve the code subspace 
are associated with the code's normalizer group, the 
group of phase space translations that commute with the 
code stabilizer. The generator matrix of the normalizer 
is a matrix M that can be chosen to be 

M x = A~ X M , (31) 

so that 

M x wM T = I ; (32) 



lattice, both A and (A x ) are integral matrices; there- 
fore D = D^ 1 and the standard form of A is 

A -(X »)="■ (38) 

Hence the generator matrix of a symplectically self-dual 
lattice can be chosen to be a real symplectic matrix: 
MloM t = uj. 

If the lattice is rotated, then the generator matrix is 
transformed as 

M -» MO , (39) 

where O is an orthogonal matrix. Therefore, it is conve- 
nient to characterize a lattice with its Gram matrix 

G = MM T , (40) 

which is symmetric, positive, and rotationally invariant. 
In the case of a symplectically self-dual lattice, the Gram 
matrix G can be chosen to be symplectic, and two sym- 
plectic Gram matrices G and G' describe the same lattice 
if 

G' = RGR T , (41) 



and 

(M x ) u (M x f = (A- l f . (33) 

We will refer to the lattice generated by M as the sym- 
plectic dual C 1 - of the lattice C. 

Another matrix that generates the same lattice as M 
(and therefore defines a different set of generators for the 
same stabilizer group) is 

M' = RM , (34) 

where R is an integral matrix with deti? = ±1. This 
replacement changes the matrix A according to 

A -> RAR T . (35) 

By Gaussian elimination, an R can be constructed such 
that 

A =( D o) • ( 36 ) 

and 

where D is a positive diagonal integral N x N matrix. 
In the important special case of a symplectically self- dual 



where R is symplectic and integral. Therefore, the mod- 
uli space of symplectically self-dual lattices in 2N dimen- 
sions can be represented as 

A N = H(2N)/Sp(2N, Z) , (42) 

where H{2N) denotes the space of real symplectic posi- 
tive 2N x 2N matrices of determinant 1. The space An 
can also be identified as the moduli space of principally 
polarized abelian varieties in complex dimension ./V [ [I7[ . 

The encoded operations that preserve the code space 
but act trivially within the code space comprise the quo- 
tient group £ x /£. The order of this group, the ratio of 
the volume of the unit cell of C to that of C^, is m 2 , 
where m is the dimension of the code space. The vol- 
ume of the unit cell of C is | det M\ = | det A] 1 / 2 and the 
volume of the unit cell of C 1 - is | det M ± \ = \ det A^ 1 / 2 ; 
therefore the dimension of the code space is 

m = |Pf A\ = |detM| =detD , (43) 

where Pf A denotes the Pfafnan of A, the square root of 
its determinant. Thus, a symplectically self-dual lattice, 
for which | det M\ = | det M ± \ = 1, corresponds to a code 
with a one-dimensional code space. Given a 27V x 27V 
generator matrix M of a symplectically self-dual lattice, 
we can rescale it as 

M -> VXM , (44) 
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where A is an integer, to obtain the generator matrix of 
a symplectically integral lattice corresponding to a code 
of dimension 



m = X N . 

The rate of this code, then, is 
R = log 2 A . 



(45) 



(46) 



When an encoded state is subjected to the Gaussian 
quantum channel, a phase space displacement 



(47) 



is applied. To diagnose and correct this error, the eigen- 
values of all stabilizer generators are measured, which 
determines the value of (£ 9 ,£ P ) modulo the normalizer 
lattice C . To recover, a displacement of minimal length 
is applied that returns the stabilizer eigenvalues to their 
standard values, and so restores the quantum state to the 
code space. We can associate with the origin of the nor- 
malizer lattice its Voronoi cell, the set of points in R 2N 
that are closer to the origin than to any other lattice site. 
Recovery is successful if the applied displacement lies in 
this Voronoi cell. Thus, we can estimate the likelihood of 
a decoding error by calculating the probability that the 
displacement lies outside the Voronoi cell. 



IV. ACHIEVABLE RATES FROM EFFICIENT 
SPHERE PACKINGS 



One way to establish an achievable rate for the Gaus- 
sian quantum channel is to choose a normalizer lattice 
C whose shortest nonzero vector is sufficiently large. In 
this Section, we calculate an achievable rate by demand- 
ing that the Voronoi cell surrounding the origin contain 
all typical displacements of the origin in the limit of large 
N. In Sec. V, we will use a more clever argument to im- 
prove our estimate of the rate. 

The volume of a sphere with unit radius in n dimen- 
sions is 



i/2 



and from the Stirling approximation we find that 



V n < 



2ire 



i/ 2 



(49) 



It was shown by Minkowski JTf| that lattice sphere pack- 
ings exist in n dimensions that fill a fraction at least 



l/2( n ~ 1 * > of space. Correspondingly, if the lattice is cho- 
sen to be unimodular, so that its unit cell has unit vol- 
ume, then kissing spheres centered at the lattice sites can 
be chosen to have a radius r„ such that 



or 



V n (r„) n > 2-("" 1 ) 



(50) 



(51) 



This lower bound on the efficiency of sphere packings 
has never been improved in the nearly 100 years since 
Minkowski's result. More recently, Buser and Sarnak |l7| ] 
have shown that this same lower bound applies to lattices 
that are symplectically self-dual. 

Now consider the case of n — 2V-dimensional phase 
space. For sufficiently large n, the channel will apply 
a phase space translation by a distance which with high 
probability will be less than \J n(a 2 + e), for any positive 
e. Therefore, a code that can correct a shift this large 
will correct all likely errors. What rate can such a code 
attain? If the code is a lattice stabilizer code, and the 
dimension of the code space is m, then the unit cell of 
the code's normalizer lattice has volume 



A = (2-kH) 

m 



N 



(52) 



Nonoverlapping spheres centered at the sites of the 
normalizer lattice can be chosen to have radius r = 
\Jn{o 2 + e), where 

^/2 



/27TC 

\ n 



(n(a 2 + e)) n ' 2 > - ■ 2"" ■ {2^h) n ' 2 , (53) 
m 



m > 



4e(cr 2 +e) 



N 



(54) 



The error probability becomes arbitrarily small for large 
N if eq. (|I|) is satisfied, for any positive e. We conclude 
that the rate 



R 



1 

N 



log 2 m = log 2 



4e<7' 



(55) 



(48) is achievable, provided h/Aea 2 > 1. However, as noted 
in Sec. Ill, the rates that can be attained by this con- 
struction (rescaling of a symplectically self-dual lattice) 
are always of the form log 2 A, where A is an integer. 



V. IMPROVING THE RATE 



The achievable rate found in eq. 
short of the coherent information eq. 



falls two qubits 
(|23J). We will now 
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show that this gap can be closed by using tighter esti- 
mates of the error probability. We established eq. (|55| ) 
by filling phase space with nonoverlapping spheres, which 
is overly conservative. It is acceptable for the spheres to 
overlap, as long as the overlaps occupy an asymptoti- 
cally negligible fraction of the total volume, as suggested 
in Fig. 0. 

Our improved estimate applies another result obtained 
by Buser and Sarnak |l7]. They note that the moduli 
space of symplectically self-dual lattices is compact and 
equipped with a natural invariant measure. Therefore, it 
makes sense to consider averaging over all lattices. De- 
note by (•) the average over all symplectically self-dual 
lattices with specified dimension n — 2N, and let f(x) de- 
note an integrable rotationally-invariant function of the 
vector x (that is a function of the length |x| of x). Then 
Buser and Sarnak [|~7| show that 

( E /(*)) = / /(*) d " x ■ ( 56 ) 

x£C\{0} J 

(Note that the sum is over all nonzero vectors in the 
lattice C.) It follows that there must exist a particular 
symplectically self-dual lattice C such that 




The statement that a unimodular lattice exists that sat- 
isfies eq. ([37|) is the well-known Minkowski-Hlawka the- 
orem Buser and Sarnak established the stronger 
result that the lattice can be chosen to be symplectically 
self-dual. 




(a) (b) 



FIG. 1. Two ways to estimate the rate achieved by a lattice 
code. Each site of the normalizer lattice has a Voronoi cell 
(represented here by a square) containing all points that are 
closer to that site than any other site. Displacements that 
move a site to a position within its Voronoi cell can be cor- 
rected. The volume of the Voronoi cell determines the rate 
of the code. In (a), the ball containing typical displacements 
lies within the cell, so that the error probability is small. In 
(6), the ball of typical displacements is not completely con- 
tained within the cell, but the region where neighboring balls 
overlap (shown in black) has a small volume, so that the error 
probability is still small. 



We can use this result to bound the probability of a de- 
coding error, and establish that a specified rate is achiev- 
able. Our argument will closely follow de Buda [ p0[ , 
who performed a similar analysis of lattice codes for the 
Gaussian classical channel. However, the quantum case 
is considerably easier to analyze, because we can avoid 
complications arising from the power constraint ]2l|-|23|| . 

A decoding error occurs if the channel displaces the 
origin to a point outside the Voronoi cell centered at the 
origin. The Voronoi cell has a complicated geometry, so 
that the error probability is not easy to analyze. But 
we can simplify the analysis with a trick pp| ]. Imagine 
drawing a sphere with radius 

a = vMo- 2 + e) (58) 

around each lattice site, where e > 0; this value of a is 
chosen so that the typical displacement introduced by the 
channel has length less than a; the probability of a shift 
larger than a thus becomes negligible for large n. It may 
be that these spheres overlap. However, a vector that is 
contained in the sphere centered at the origin, and is not 
contained in the sphere centered at any other lattice site, 
must be closer to the origin than any other lattice site. 
Therefore, the vector is contained in the origin's Voronoi 
cell, and is a shift that can be corrected successfully. (See 
Fig. 1.) 

Hence (ignoring the possibility of an atypical shift by 
£ > a) we can upper bound the probability of error by 
estimating the probability that the shift moves any other 
lattice site into the sphere of radius a around the origin. 
We then find 

Perror < ^ / P ( X ~ T )^ T ' ( 59 ) 

xec^\{o} J \ r ^ a 
where P(£) denotes the probability of a displacement by 

The Buser-Sarnak theorem |l7|] tells us that there ex- 
ists a lattice whose unit cell has volume A, and which is 
related by rescaling to a symplectically self-dual lattice, 
such that 

Perror < ~ J <Tx J ^ P( X ~ r)d n r | (60) 

by interchanging the order of integration, we find that 

Perror < ^ • V n • O n , (61) 

the ratio of the volume of the n-dimensional sphere of 
radius a to the volume of the unit cell. 

Now the volume A of the unit cell of the normalizer 
lattice C , and the dimension m of the code space, are 
related by 
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A = {2-Kh) N mT 1 = (27rft-2" 



(62) 



where i? is the rate, and we may estimate the volume of 
the sphere as 



K-a"<(^)" /2 (^ + e))" /2 



where n — 2N. Thus we conclude that 



iV 



(63) 



(64) 



Therefore, the error probability becomes small for large 
TV for any rate R such that 



R < log 2 



e(tr 2 + e) 



(65) 



where e may be arbitrarily small. We conclude that the 
rate 



R = lot 



2 ? 



(66) 



if a sufficiently small fraction of the N qudits are dam- 
aged. Let us compare the rates achieved by concatenated 
codes to the rates achieved with codes derived from effi- 
cient sphere packings. 

We analyze the effectiveness of concatenated codes in 
two stages. First we consider how likely each of the N 
qudits is to sustain damage if the underlying oscillator is 
subjected to the Gaussian quantum channel. The area 
of the unit cell of the two-dimensional square normalizer 
lattice that represents the encoded operations acting on 
the qudit is 2nh/d, and the minimum distance between 
lattice sites is S — ^2irh/d. A displacement of q by a ■ <5, 
where a is an integer, is the operation X a acting on the 
code space, and a displacement of p by b-S is the operation 
Z b , where A and Z are the Pauli operators acting on the 
qudit; these act on a basis {\j), j = 0, 1, 2, . . . , d— 1} for 
the qudit according to 



X 

z 



\j + 1 (mod d)) 
w 3 \j) , 



(69) 



where u> = exp(27ri / d) . 

Shifts in p or q can be corrected successfully provided 
that they satisfy 



is achievable in the limit N — > oo, provided that h/ea 2 > 
1. This rate matches the optimal value eq. (^) of the 
one-shot coherent information for Gaussian inputs. We 
note, again, that the rates that we can obtain from rescal- 
ing a symplectically self-dual lattice are restricted to 
R = log 2 A, where A is an integer. Thus for specified 
a 2 , the achievable rate that we have established is really 
the maximal value of 



R = log 2 A , A 6 Z , 
such that the positive integer A satisfies 

n 



A < 



(67) 



(68) 



eo" 1 



VI. ACHIEVABLE RATES FROM 
CONCATENATED CODES 



Another method for establishing achievable rates over 
the Gaussian quantum channel was described in p4[ , 
based on concatenated coding. In each of N "oscilla- 
tors" described by canonical variables pi and a d- 
dimensional system ( "qudit" ) is encoded that is pro- 
tected against sufficiently small shifts in pi and qi. The 
encoded qudit is associated with a square lattice in 2- 
dimensional phase space. Then a stabilizer code is con- 
structed that embeds a fc-qudit code space in the Hilbert 
space of N qudits; these k encoded qudits are protected 



\Aq\ < 5/2 



|Ap| < 6/2 = 



(70) 



If the shifts in q and p are Gaussian random variables 
with variance a 2 , then the probability that a shift causes 
an uncorrectable error is no larger than the probability 
that the shift exceeds WTfh/2d, or 



Px,Pz < 



V2 



7TO"' 



dxe- x2 ' 2 ° 



eric^nh/Ada 2 ) , 



(71) 



where erfc denotes the complementary error function. 
Here px is the probability of an "A error" acting on 
the qudit, of the form X a for a ^ (mod d), and pz 
denotes the probability of a 11 Z error" of the form Z for 
b ^ (mod d). The A and Z errors are uncorrelated, and 
errors with a, b = ±1 are much more likely than errors 
with |a|, \ b\ > 1. By choosing d ~ ti/o- 2 , we can achieve 
a small error probability for each oscillator. 

The second stage of the argument is to determine the 
rate that can be achieved by a qudit code if px , Pz satisfy 
eq. (|7l|) . We will consider codes of the Calderbank-Shor- 
Steane (CSS) type, for which the correction of A errors 
and Z errors can be considered separately |2^j2^]. A 
CSS code is a stabilizer code, in which each stabilizer 
generator is either a tensor product of Ps and powers of 
Z (measuring these generators diagnoses the A errors) or 
a tensor product of 7's and powers of A (for diagnosing 
the Z errors). 
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We can establish an achievable rate by averaging the 
error probability over CSS codes; we give only an infor- 
mal sketch of the argument. Suppose that we fix the 
block size N and the number of encoded qudits k. Now 
select the generators of the code's stabilizer group at ran- 
dom. About half of the N — k generators are of the Z 
type and about half are of the X type. Thus the number 
of possible values for the eigenvalues of the generators of 
each type is about 



rf 4(JV-fe) 



(72) 



Now we can analyze the probability that an uncorrectable 
X error afflicts the encoded quantum state (the proba- 
bility of an uncorrectable Z error is analyzed in exactly 
the same way). Suppose that X errors act independently 
on the N qudits in the block, with a probability of error 
per qudit of px- Thus for large N, the typical number of 
damaged qudits is close to px ■ N. A damaged qudit can 
be damaged in any of d — 1 different ways (X a , where 
a = 1, 2, . . . , (d — 1)). We will suppose, pessimistically, 
that all d — 1 shifts of the qudit are equally likely. The 
actual situation that arises in our concatenated coding 
scheme is more favorable - small values of \a\ are more 
likely - but our argument will not exploit this feature. 

Thus, with high probability, the error that afflicts the 
block will belong to a typical set of errors that contains 
a number of elements close to 



typ 



( N 

\Np x 



Vrf - 1)** ~ ( lN(H d (p x )+p x log d (d-l)) ^ 



where 



H d (p) = -p\og d p- (1 - p) log d (l - p) 



(73) 



(74) 



If a particular typical error occurs, then recovery will suc- 
ceed as long as there is no other typical error that gener- 
ates the same error syndrome. It will be highly unlikely 
that another typical error has the same syndrome as the 
actual error, ^provided that the number of possible error 
syndromes d^ N ~ k > is large compared to the number of 
typical errors. Therefore, the X errors can be corrected 
with high probability for 



1/ * 
2 V N 



> 



1 

N 



\og d N t 



typ 



H d {p x )+Pxlog d (d-l) , (75) 



or for a rate R d in qudits satisfying 

Rd = jj < 1 - 2H d ( Px ) - 2 PX log d (d - 1) (76) 

Similarly, the Z errors can be corrected with high prob- 
ability by a random CSS code if the rate satisfies 



R d < 1 - 2H d ( Pz ) - 2p z \og d {d - 1) 

Converted to qubits, the rate becomes 
R = log 2 d ■ R d 



(77) 



(78) 



Under these conditions, the probability of error averaged 
over CSS codes becomes arbitrarily small for N large. It 
follows that there is a particular sequence of CSS codes 
with rate approaching eq. ( f76f|78| ), and error probability 
going to zero in the limit N — > oo. 

For given a 2 , the optimal rate that can be attained by 
concatenating a code that encodes a qudit in a single os- 
cillator with a random CSS code, is found by estimating 
px and pz using eq. ( |7l| ) and then choosing d to max- 
imize the rate R given by eq. (f7(||7|). The results are 
shown in Fig. |^. This rate (in qubits) can be expressed 
as 



R = log 2 (C 2 h/a 2 



(79) 



where C 2 is a slowly varying function of a 2 /H plotted in 
Fig. H It turns out that this rate is actually fairly close 
to log 2 d; that is, the optimal dimension d of the qu- 
dit encoded in each oscillator is approximately C 2 h/a 2 . 
With this choice for d, the error rate for each oscillator is 
reasonably small, and the random CSS code reduces the 
error probability for the encoded state to a value expo- 
nentially small in N at a modest cost in rate. The rate 
achieved by concatenating coding lies strictly below the 
coherent information Iq, but comes within one qubit of 
Iq for cr 2 > 1.88 x 1(T 4 . 




FIG. 2. Rates achieved by concatenated codes, compared 
to the one-shot coherent information optimized over Gaus- 
sian input states. Here a is the standard deviation of the 
magnitude of the phase-space displacement introduced by the 
channel, in units with h = 1. 

Both the concatenated codes and the codes de- 
rived from efficient sphere packings are stabilizer codes, 
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and therefore both are associated with lattices in 2N- 
dimensional phase space. But while the sphere-packing 
codes have been chosen so that the shortest nonzero vec- 
tor on the lattice is large relative to the size of the unit 
cell, the concatenated codes correspond to sphere pack- 
ings of poor quality. For the concatenated codes, the 
shortest vector of the normalizer lattice has length £, 
where 



I 2 = 2irh/d 



(80) 



and the rate R is close to log 2 d. The efficient sphere 



packings have radius r = 1/2 close to 



8Nh 



(81) 



Hence, if we compare sphere-packing codes and concate- 
nated codes with comparable rates, the sphere-packing 
codes have minimum distance that is larger by a factor 
of about y/4N/Tre. The concatenated codes achieve a 
high rate not because the minimum distance of the lat- 
tice is large, but rather because the decoding procedure 
exploits the hierarchical structure of the code. 



0.4 n 

0.3- 
0.2 



0.1- 



.001 



.01 



FIG. 3. The slowly varying function C , defined by 
R — log 2 (C 2 /<T 2 ), where R is the rate achievable with con- 
catenated codes. Units have been chosen such that h = 1. 
The horizontal lines are at C 2 = 1/e, corresponding to a rate 
equal to the coherent information, and at C 2 = l/2e, corre- 
sponding to one qubit below the coherent information. 



VII. THE CLASSICAL GAUSSIAN CHANNEL 

We have found that quantum stabilizer codes based on 
efficient sphere packings can achieve rates for the Gaus- 
sian quantum channel that match the one-shot coherent 
information, and that concatenated codes achieve rates 
that are below, but close to, the coherent information. 
Now, as an aside, we will discuss the corresponding state- 
ments for the classical Gaussian channel. We will see, in 



particular, that concatenated codes achieve rates that are 
close to the classical channel capacity. 

Shannon's expression for the capacity of the classical 
Gaussian channel can be understood heuristically as fol- 
lows If the input signals have average power P, 
which is inflated by the Gaussian noise to P + a 2 , then 
if n real variables are transmitted, the total volume oc- 
cupied by the space of output signals is the volume of a 
sphere of radius y/n(P + a 2 ), or 



total volume = V n ■ (n(P + a 2 )) 



n/2 



(82) 



We will decode a received message as the signal state 
that is the minimal distance away. Consider averaging 
over all codes that satisfy the power constraint and have 
to signals. When a message is received, the signal that 
was sent will typically occupy a decoding sphere of radius 
■J (n(<j 2 + e) centered at the received message, which has 
volume 



decoding sphere volume = V n ■ (n(a 2 + e)) 



i/2 



(83) 



A decoding error can arise if another one of the to signals, 
aside from the one that was sent, is also contained in 
the decoding sphere. The probability that a randomly 
selected signal inside the sphere of radius y/n(P + cr 2 ) 
is contained in a particular decoding sphere of radius 
\J n(a 2 + e) is the ratio of the volume of the spheres, so 
the probability of a decoding error can be upper bounded 
by m times that ratio, or 



< to ■ 



P 



i/2 



>2H 



P 



i/2 



(84) 



where R is the rate of the code. If the probability of er- 
ror averaged over codes and signals satisfies this bound, 
there is a particular code that satisfies the bound when 
we average only over signals. If P C rror < 8 when we av- 
erage over signals, then we can discard at most half of 
all the signals (reducing the rate by at most l/n bits) to 
obtain a new code with P e rror < 25 for all signals. Since 
e can be chosen arbitrarily small for sufficiently large n, 
we conclude that there exist codes with arbitrarily small 
probability of error and rate R arbitrarily close to 



1 / P 
C=-lo g2 



(85) 



which is the Shannon capacity. Conversely, for any rate 
exceeding C, the decoding spheres inevitably have non- 
negligible overlaps, and the error rate cannot be arbitrar- 
ily small. 

Suppose that, instead of Shannon's random coding, we 
use a lattice code based on an efficient packing of spheres. 
In this case, the power constraint can be imposed by in- 
cluding as signals all lattice sites that are contained in an 
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n-dimensional ball of radius ynP : and the typical shifts 
by distance V no 2 must be correctable. Thus decoding 

are to be packed into a sphere 



spheres of radius v 
of total radius \J n(P + a 2 ). Suppose that the lattice is 
chosen so that nonoverlapping spheres centered at the 
lattice sites fill a fraction at least 2 _ (™ _1 ) of the total 
volume; the existence of such a lattice is established by 
Minkowski's estimate |18 . Then the number m of signals 
satisfies 



m-V n - (na 2 ) n/2 > 2-(™~ 1 ) • V„ ■ (n(P + a 2 )) 



i/2 



or 



m > T 



n/2 



(86) 



(87) 



corresponding to the rate 

~ 1 , 1 , ( P 

R = - ■ log 2 m = - log 2 1 + — 

n 2 V cr 



which is one bit less than the Shannon capacity. 

Much as in the discussion of quantum lattice codes 
in Sec. [v| an improved estimate of the achievable rate 
is obtained if we allow the decoding spheres to overlap 
p(j| p3[ . In fact, there are classical lattice codes with rate 
arbitrarily close to the capacity, such that the probability 
of error, averaged over signals, is arbitrarily small p3[ . 
Unfortunately, though, because of the power constraint, 
the error probability depends on which signal is sent, and 
the trick of deleting the worst half of the signals would 
destroy the structure of the lattice. Alternatively, it can 
be shown that for any rate 



R<^log 2 (P/a 2 ) , 



(89) 



there are lattice codes with maximal probability of er- 
ror that is arbitrarily small p0| . This achievable rate 
approaches the capacity for large P/cr 2 . 

Now consider the rates that can be achieved for the 
Gaussian classical channel with concatenated coding. A 
d-state system (dit) is encoded in each of n real variables. 
If each real variable takes one of d possible values, with 
spacing 2Ax between the signals, then a shift by Aa; can 
be corrected. By replacing the sum over d values by an 
integral, which can be justified for large d, we find an 
average power per signal 



P 



1 



2dAx 



dAx 



dAx 



x z dx = ]-{dAxf 



(90) 



thus the largest correctable shift can be expressed in 
terms of the average power as 



For the Gaussian channel with mean zero and variance 
a 2 , the probability p of an error in each real variable 
transmitted is no larger than the probability of a shift by 
a distance exceeding Ax, or 



p < erfc (v/3P/2d 2 cr : 



(92) 



where erfc denotes the complementary error function. 

We reduce the error probability further by enco ding 
k < n dits in the block of n dits. Arguing as in Sec. VI, 
we see that a random code for dits achieves an asymptotic 
rate in bits given by 



R 



\og 2 d 



(l-H d (p)~plog d (d-l)) 



(93) 



Given cr 2 , using the expression eq. (£21) for p, and choos- 
ing d to optimize the rate in eq. (|93[), we obtain a rate 
close to the Shannon capacity, as shown in Fig. [|. As 
for the concatenated quantum code, the rate of the con- 
1 ; (88) catenated classical code is close to log 2 d, where d ~ 

C(cr 2 ) • yj P/a 2 , and C(cr 2 ) is a slowly varying function. 




Aa; = V3P/d 



(91) 



FIG. 4. Rates for the Gaussian classical channel achievable 
with concatenated codes, compared to the Shannon capacity. 
Here a is the standard deviation of the displacement, in units 
with the power P = 1. 



VIII. CONCLUSIONS 



We have described quantum stabilizer codes, based 
on symplectically integral lattices in phase space, that 
protect quantum information carried by systems de- 
scribed by continuous quantum variables. With these 
codes, we can establish lower bounds on the capacities of 
continuous- variable quantum channels. 

For the Gaussian quantum channel, the best rate we 
know how to achieve with stabilizer coding matches the 
one-shot coherent information optimized over Gaussian 
inputs, at least when the value of the coherent infor- 
mation is log 2 of an integer. That our achievable rate 
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matches the coherent information only for isolated val- 
ues of the noise variance a 2 seems to be an artifact of 
our method of analysis, rather than indicative of any in- 
trinsic property of the channel. Hence it is tempting to 
speculate that this optimal one-shot coherent informa- 
tion actually is the quantum capacity of the channel. 

Conceivably, better rates can be achieved with nonad- 
ditive quantum codes that cannot be described in terms 
of symplectically integral lattices. We don't know much 
about how to construct these codes, or about their prop- 
erties. 

In the case of the depolarizing channel acting on qubits, 
Shor and Smolin discovered that rates exceeding the 
one-shot coherent information could be achieved. Their 
construction used concatenated codes, where the "outer 
code" is a random stabilizer code, and the "inner code" 
is a degenerate code with a small block size El|. The 
analogous procedure for the Gaussian channel would be 
to concatenate an outer code based on a symplectically 
integral lattice with an inner code that encodes one logi- 
cal oscillator in a block of several oscillators. This inner 
code, then, embeds an infinite-dimensional code space 
in a larger infinite-dimensional space, as do codes con- 
structed by Braunstein |^6j and Lloyd and Slotine [ p7[ . 
However, we have not been able to find concatenated 
codes of this type that achieve rates exceeding the one- 
shot coherent information of the Gaussian channel. 
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